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Electronic spin-orbit coupling (SOC) is essential for various newly discovered phenomena in 
condensed-matter systems. In particular, one-dimensional topological heterostructures with SOC 
have been widely investigated in both theory and experiment for their distinct transport signa¬ 
tures indicating the presence of emergent Majorana fermions. However, a general framework for 
the SOC-affected transport in superconducting heterostructures, especially with the consideration 
of interfacial effects, has not been developed even regardless of the topological aspects. We hereby 
provide one for an effectively one-dimensional superconductor-normal heterostructure with nonuni¬ 
form magnitude and direction of both Rashba and Dresselhaus SOC as well as a spin-orbit-active 
interface. We extend the Blonder-Tinkham-Klapwijk treatment to analyze the current-voltage re¬ 
lation and obtain a rich range of transport behaviors. Our work provides a basis for characterizing 
fundamental physics arising from spin-orbit interactions in heterostructures and its implications for 
topological systems, spintronic applications, and a whole variety of experimental setups. 

PACS numbers: 74.45.+c, 74.25.F-, 71.70.Ej, 74.25.fc 

I. INTRODUCTION 

The interplay between electronic spin and orbital de¬ 
grees of freedom, or spin-orbit coupling (SOC), has 
played a crucial role in various aspects of condensed- 
matter physics, including the study on semiconductors, 
ferromagnets, superconductors and materials with exotic 
orders [1-14] as well as the application on building quan¬ 
tum electronic or spintronic devices [15, 16]. Due to the 
ubiquity of SOC and superconductivity, it is of funda¬ 
mental and technological interest to understand their in¬ 
terplay, in particular, in the context of superconducting 
heterostructures. Furthermore, great interest has been 
stimulated recently in exploring topological states of mat¬ 
ter [17, 18], whose experimental realization and detec¬ 
tion strongly rely on transport or scanning tunneling mi¬ 
croscopy (STM) measurements of artificially engineered 
heterostructures with SOC. Among them, those with 
building blocks such as a semiconductor/topological- 
insulator wire with SOC and conventional superconduc¬ 
tor [19-23] (used to induce a proximity gap) or one¬ 
dimensional (ID) ferromagnet combined with a bulk su¬ 
perconductors with SOC [24-28] are of particular inter¬ 
est to explore emerging Majorana fermions [29-35]. Most 
theoretical studies have focused on either effective mod¬ 
els or employed SOC as a uniform model parameter [36- 
47]. However, the set-up of interest usually has only a 
segment or end point of the ID element in contact with 
a (bulk) superconductor due to either the experimental 
constraints or an explicit interest in studying (topologi¬ 
cal) heterostructures [23, 40, 48, 49]. 

In this paper, we develop a framework to study an (ef¬ 
fectively) ID normal- (N-) superconductor (S) setup (see 
Fig. 1) that is general enough to capture the effects of 
SOC in a wide range of (emergent) heterostructures. We 



FIG. 1. (Color online) Illustration of an effective ID N-S 
system with SOC and a spin-orbit-active interface (I). The 
vectors denote directions of Rashba aN(S) ff i/ and Dresselhaus 
i$n( s)°4: on the N (S) side, as well as interfacial Rashba Z a a y , 
Dresselhaus Z a a x , and Zeeman Zhfi-ff components. In gen¬ 
eral, the spin basis in each region can be different. Insets: 
schematic energy spectra with two isospin branches (solid and 
dashed curves, respectively) on the corresponding sides. The 
solid circles denote possible outgoing states corresponding to 
an incoming electron (empty circle) from the N side. 

allow for both Rashba [50] and Dresselhaus [51] SOC with 
arbitrary relative strengths on each side as well as at the 
interface (I), which represents the boundary region [52]. 
Furthermore, we allow the SOC to have a different mag¬ 
nitude as well as direction in N, S, and I regions. Nonuni¬ 
form SOC is in fact bound to arise for a variety of reasons 
including difference in the chemical composition, crystal¬ 
lographic direction, or direction of electric field [50, 51] 
and may be relevant also for practical implementations of 
topological quantum computing schemes [53, 54]. SOC 
in the interface can even simply arise due to a lack of 
“left-right” symmetry [55-57]. 

We start by analyzing the appropriate Bogoliubov- 
de Gennes (BdG) Hamiltonian and develop a general¬ 
ized Blonder-Tinkham-Klapwijk (BTK) [58] formalism 
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to calculate the current-voltage characteristics for our 
general setup. We phenomenologically model what we 
call the spin-orbit-active interface as a generalization of 
a spin-active interface [59-74], which itself is a gener¬ 
alization of a single Z parameter used in the original 
BTK approach. We find that the rich physics emanat¬ 
ing from non-uniform SOC and/or spin-orbit-active in¬ 
terface is not only of fundamental interest but also highly 
relevent for interpreting the results of experiments. For 
example, it can effectively lower the interfacial barrier, 
drastically change the zero-bias conductance, and deter¬ 
mine between reentrant or monotonic temperature de¬ 
pendence. Our results also underline the danger in using 
a single /^-parameter fit to extract the interface proper¬ 
ties in experiments as is standard practice. Our analysis 
offers a more accurate picture even for traditional N-S 
heterostructures or for STM setups in which a spin-orbit- 
active interface or the nonzero SOC in the materials may 
play a role that is not usually taken into account (al¬ 
though the importance of nonzero SOC in the context of 
pair-breaking effects [75-80] is appreciated). 

II. MODEL AND HAMILTONIAN 

As illustrated in Fig. 1, our effectively ID system com¬ 
prises N (x < 0) and S (x > 0) regions with an ar¬ 
bitrary SOC on each side, intersecting at I (x = 0). 
The Hamiltonian _Hj = n,s i n corresponding basis ’Fj = 

V’H i’t -V’t,. ) reads 

Hn = T z <g> (E p t + a N pa y + f3 N pcr x ), (1) 

H s =t z <8> (E p t + a s pcTy + /3 s pa x ) - A t x <8> 1, (2) 
2 

where E p = — p is the free spectrum with chemical 

potential p, ay (//,) are Rashba (Dresselhaus) SOC on 
the j side, A (taken real for convenience) is the s-wave 
(possibly proximity-induced) superconducting gap, a and 
r are Pauli matrices spanned in spin and particle-hole 
basis, respectively, and 1 is the 2x2 identity matrix. The 
spin basis on both sides can differ by an Euler rotation 


T n = /I (§) e~ iazTlNS e~ iayr1 ^e~ ilTzri ^j % 

= Uisss'&s- ( 3 ) 

The two SOCs can be combined as a complex factor 

7 =/3j + icij. (4) 

We rotate the Hamiltonians into an isospin basis Ty un¬ 
der transformation Hy = UjyHjUjji where 

1 / 1 \ 

-1 ( 5 ) 

and obtain a diagonal H N > and block-diagonal H$i as 

#N' =t z ® ( E p t - 7nP<Tz) , (6) 

Hs> = t z <8> (Ept - 7SP07) - At x (g> 1. (7) 


The spectra of Hn' and ifs' show particle and hole bands, 
denoted by r = ±. Each band has two isospin branches 
denoted by a = ±. At a given energy E, there are eight 
corresponding eigen wave functions for H^i and the other 
eight for Hs>, denoted as 


e^ ±v - +a ^ x and ^(±*£+*78)®, 

(8) 

respectively, where 


p T a = y/2m (p' — mdy 2 + tE), 

(9) 

= 2m ^ p' + mS"/ 2 + t\JE 2 — A 2 j , 

(10) 

with 


p! = p + m(7 N+7s) an( J _ 7s -7n _ 

(11) 


The particle-hole band and isospin states have represen¬ 
tations as Xn ,<T = ( ^r+dcr+ St+ 5,7- S T -S a + S T -S a - ) 
and Xs' a = u ( S T+ S a+ 8 T+ 8 a - 8 T _5 a+ 8 T S a - ) T + 
v( S-T-Sa + 8 X — 8U — 8 T j r 8 CT j r 8 t +8<j _ ) T , where u 2 = 1 — 
v 2 = 1 ^1 + ^ and 8 is the delta function. Below 

we drop H in all equations for convenience and take p' 
and i/2m/7 to be natural energy and momentum units, 
respectively (so 7 g + 7 ^ causes no qualitative change but 
energy rescaling). Note that we consider energy range in 
which the hole excitations have real momenta [81]. 

In a BTK treatment, the interface is phenomeno¬ 
logically modeled by Hamiltonian H\8(x ), which gen¬ 
erally has the same 4x4 representation as the bulk 
Hamiltonian. The matrix components reflect the in¬ 
terfacial properties as well as the physical discontinu¬ 
ity between both sides and specify transport processes 
through the interface. The original BTK model adopts 
one parameter Z 0 t z (8) 1 to describe barrier effects from 
an oxide layer or local disorder [58]. For ferromagnet- 
superconductor heterostructures, Hi can be modeled as 
a Zeeman form 1 ( 8 > Zh (n ■ a) with unit direction n = 
(sin(9cos</, sin 0 sin ^ 1 , cos 0 ), which accounts for various 
spin-active processes such as spin-flip scattering, spin- 
dependent phase shift, and spin-related Andreev reflec¬ 
tion [61-74]. In our case, we expect that the discontinuity 
of the bulk SOC and the interfacial SOC itself would play 
an active role in the transport. This effect is modeled by 
Rashba Z a T z (g> a y and Dresselhaus Zqt z (g> <j x compo¬ 
nents. To capture the most general interplay with spins, 
we incorporate all the factors above and write down our 
spin-orbit-active interface, 

H\ = t z O (Z 0 1 + Z a cr y + Zpcr x ) + 1 (8) Z h (n ■ a) (12) 

The two SOCs can be combined as Zp + iZ a = Z 7 e l,pi . 
In general, the spin basis of Hi and H n can differ by 
another Euler rotation 

Um = (l 8> e~ icr ^ . (13) 
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We write the rotated H\ in a suggestive form as 


/ 

Zq — Z\ 

Z 2 e~ i < 

0 

0 

) 



Z 2 e* 

Z 0 + Zi 

0 

0 

f° + { 


0 

0 

. CO 

N 

-Z 4 e- i(c+c,) 


L \ 

V 

0 

0 

-Z 4 e^ + ^ 

—Z 0 + Z 3 

J 


(14) 


while the second one can be obtained by integrating the 
BdG equation over an infinitesimal interval across the 
interface, 


1, 


+U m H 1 U^ 1 S{x)\^dx = 0, 


(19) 


There are relations Z\ + Z\ + Z\ + Zf = Z\ + Z^ and 
Z\ + Z\ — Z\ — Z\ = 2Z^Zh sind cos (</> — <pi) that are 
independent of the Euler angles. Comparing Eq. (12) 
and Eq. (14), one can see that the three Euler angles 
effectively generate one additional variable. 

In brief, the key ingredients on our system, the bulk 
SOCs ctjCTy and /3jcr x as well as the interfacial parameters 
Z a (j y , Zpa Xl and Zhfi-a , can be characterized by differ¬ 
ent vectors in spin space in the corresponding regions, 
as illustrated in Fig. 1. Below we study the system’s 
transport properties as a function of these vectors. 

III. BTK CALCULATIONS 


where 6{x) is the step function and \H((P) = Tr/i/O). 
We carefully keep the Hamiltonian Hermitian when it is 
expressed using the step function [83] . 

The probability current J is calculated from its def¬ 
inition = —d x J. For each incoming wave 'F j r n ' T , 

we calculate J corresponding to different scattering pro¬ 
cesses normalized by the incoming current and obtain 
combined currents carried by the in- and cross-branch 
normal (Andreev) reflections together, J™ (J^), as 
well as those carried by all the transmissions together, 
Jj a (= 1 — J™ — due to the probability conserva¬ 
tion). Following the standard BTK treatment, the net 
charge current / induced by a voltage drop V across the 
junction can be evaluated as 


Here we apply the BTK treatment to compute current- 
voltage relations of the system. Considering an incoming 
wave 


,T,r,cr T.cr 

^in = X N ' e 


i(rpJx+<T7 n) 


(15) 


with energy E on the N side that propagates toward the 
interface (positive group velocity) and scatters through, 
we incorporate all possible outgoing waves and write 
down the wave functions on the N and S sides, as 


0 W7n 


= J7nn' [m + E c 
X [KXn' 17 

Tr = C^ns^ss' E e 

Z -'(7 

x (c,xr^+iX ^ e ~ ih ’ x ). 


e --~+<UXn^" x ) 

la'ys 


(16) 

(17) 


respectively (all wave functions are represented in the 
real spin basis of H^). Because reflections (transmis¬ 
sions) should have the group velocity direction opposite 
to (same as) the incoming wave, only eight states are 
considered (see insets in Fig. 1). For an incoming elec¬ 
tron (r = +), the amplitudes b, a , c, and d correspond to 
normal reflection, Andreev reflection [82], quasi-particle 
transmission and quasi-hole transmission, respectively, 
while for an incoming hole (—), a and b reverse their 
roles. The subscript a of the amplitudes describes in¬ 
branch (cross-branch) processes if its sign is the same as 
(opposite to) the incoming wave. These amplitudes can 
be determined by two boundary conditions at the inter¬ 
face. The first one is the continuity of the wavefunction 


*l(0) = 'MO), (18) 


Aer 


{(' 


+J^[f(E)-f(E + reV)]} 1 


f(E 


reV) - f(E) 

( 20 ) 


where A is a constant associated with density of states, 
Fermi velocity, as well as an effective cross-sectional area, 
and f(E) = [exp(E/kBT) + l ] _1 is the Fermi distribution 
function at temperature T. We compute the I-V relation 
and normalized differential conductance (NDC) G/Go, 
where G = ^ and Go = dI ^7°' > is a reference value 
when the S side is normal. 


IV. RESULTS 

We first find that Jt,o' / ' AR ^ T and hence I(V) are 
independent of the Euler angles Uivs• Therefore, parame¬ 
ters such as </?n and <ps that rely on the relative spin coor¬ 
dinate (RSC) between both sides should play no role on 
the transport either. We numerically confirm this inde¬ 
pendence by varying the ratio of Rashba and Dresselhaus 
SOC on each side. For the interface, I(V) is independent 
of 77 ^j if one of Zh and Z 7 in Eq. (12) is zero, so n or 
ipi has no effect given Z 1 = 0 or Zh = 0, respectively. 
We find that the scattering amplitude and probability 
current for each channel sensitively vary with RSC, but 
they compensate in the summation for I(V). This im¬ 
plies that RSC may affect I(V) as a result of interference 
in a multichannel or multiterminal system [73, 74] and 
also open interesting possibilities for spintronics. 

In Fig. 2(a) we plot zero-bias NDC in the Zq-Z 7 plane 
for a purely spin-orbit-active interface {Zh = 0). We see 
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FIG. 2. (Color online) (a) [(b)] Zero-bias NDC (scaled by 
colors in bar graph) in the Zo-Z 7 ( Zh ) plane given Zh (Z-y) = 
0. Data are for A = 0.002, S'y 2 = 0, and T = 0. 


that NDC reaches a high value when Z 1 is close to Z 3 
and exhibits a symmetry under the exchange between Zo 
and Zj. To explain this, one can Euler rotate Hi to its 
eigenspin basis. The eigenvalues can be regarded as a 
set of characteristic barrier strengths \Zq ± Z 7 \ that de¬ 
termine the BTK results. Such set is invariant under 
Z 0 Z 1 and hence leads to the symmetry. To under¬ 
stand the NDC peaks, one can look at the scattering of 
the eigenspin states. At Z 7 = Zq, electrons of one spin 
direction and holes of the opposite both see a clean inter¬ 
face |Z 0 -Z 7 | =0, which maximizes the Andreev current 
as well as NDC. In other words, the presence of interfa¬ 
cial SOC can effectively lower the original BTK barrier 
(\Z 0 — Z 7 1 < Zo) [84]. In Fig. 2(b) we plot NDC in the 
Zo~Zh plane for a purely spin-active interface Z-y = 0 for 
comparison. We see that there is no symmetry under 
Z 0 Zh and no conductance peak along the Z 0 = Zh 
line. Such differences illustrate the interfacial SOC effects 
that the original spin-active picture does not capture. 

With the coexistence of Zeeman and SOC effects at the 
interface, the Euler angles rj^ are no longer irrelevant. 
We consider the general Hamiltonian of Eq. (14) and find 
that I(V) is independent of £. The relevant variables are 
the five strength parameters Zo,i, 2 , 3,4 and the phase dif¬ 
ference (' between off-diagonal elements of the particle 
and hole blocks. Notice that the role of is special: (1) 
it comes from the interplay between Zh and Z-y and (2) 
it does not alter the eigenvalues of the interfacial Hamil¬ 
tonian. Therefore, its effects on the transport can be 
attributed to the interference between particle and hole 
channels. In Figs. 3(a)-3(c), we plot NDC vs V at various 
(/ and Zq (we set Z 3 = Z 3 = 0 and Z 2 = Z± = 0.5 for il¬ 
lustrating the salient features from the interplay between 
and Zo). At Z 0 = 0 (a), the curves are all above 1 
and show a qualitative change from center-dent to center- 
peak types as C goes through 0 (diamonds), 0.257T (tri¬ 
angles), 0.57T (squares), 0.757T (circles), and n (crosses). 
In this half period, the zero-bias NDC monotonically in¬ 
creases with . From (’ = n to = 2n, the deformation 



FIG. 3. (Color online) (a)-(c) Interfacial phase effects on 
NDC for a general spin-orbit-active interface described by 
Eq. (14) with Zo = 0, 0.66, and 1.5, respectively. Curves with 
diamonds, triangles, squares, circles, and crosses correspond 
to (' — 0, 0.257T, 0.57T, 0.75-7T, and tt, respectively. The other 
relevant parameters are set as Z\ = Z 3 = 0, Z 2 = Z± = 0.5, 
A = 0.002, <5y 2 = 0, and T = 0.15A. (d)-(f) Charge im¬ 
balance effects on NDC at cases of = 0, 0.65-7T, and 7 r, 
respectively. Curves with diamonds, triangles, squares, cir¬ 
cles, and crosses correspond to mS'y 2 = —0.975, —0.5, 0, 0.5, 
and 0.975, respectively. The other parameters are the same 
as (a)-(c) except Zo = 0 and S-y 2 varies. 


of curves completes the other half period and reverses 
back to the = 0 case. As Zo increases, the center-peak 
curves drastically deform toward the center-dent type, 
and the trend of zero-bias NDC vs also reverses. At 
Zq = 0.66 (b), the (' = n curve has the lowest zero-bias 
NDC ss 1. At Zq = 1.5 (c), the center-dent curves remain 
and are well below 1 at low bias due to the suppression of 
transmissions and Andreev reflections (tunneling limit in 
the BTK model). The role of £' is not significant in this 
case. We turn to show one way to independently tune £' 
via the tuning of parameters in Eq. (12). Assuming the 
same spin basis on the N side and the interface (r/^j 1 = 0) 
and the interfacial Zeeman components tuned as 9 = 7r/2, 
Zh cos (j) = Z-y sin 2 ipi/ cos ipi, and Zh sin <j) = —Z^siiupi, 
we obtain Z\ = Z 3 = 0, Z 2 = Z 4 = Z^/costpi, and 

= —2ipi by equating Eqs. (12) and (14). In this case, 

is associated with the ratio of Rashba and Dresselhaus 
components at the interface with properly controlled Zee- 
man components. 

All the results above are for cases of balanced charge 
carriers between N and S sides (tn = 7s)- In Figs. 3(d)- 
3(f) we explore the effects of nonuniform SOC induced 
imbalance (dy 2 ^ 0 resulting in Fermi momentum mis¬ 
match) and its interplay with the interfacial parameters. 
At C = 0, the curves are all of the center-dent type 
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FIG. 4. (Color online) (a) Zero-bias NDC vs T at various 
(,' [convention and setting are the same as Fig. 3(a)]. (b) 

Zero-bias NDC vs T at various 8 y 2 [same as Fig. 3(e)]. 


and roughly display symmetry between positive and neg¬ 
ative imbalance. The zero-bias NDC is far below 1 in the 
highly imbalanced cases. As (' increases, the Sj 2 < 0 
curves deform more drastically than the <5y 2 > 0 ones. 
At = 0.657T, the m5 ^ 2 = —0.5 curve develops a plateau 
around V = 0, indicating an incipience of a center peak. 
The m5"/ 2 = —0.975 curve develops minimums around 
eV = A. At (' = 7r, the mdj 2 = —0.5 and 0 curves 
show a center-peak structure, and the double minimums 
of the mS'y 2 = —0.975 curves become more significant. 
The range of zero-bias NDC as a function of 1 5y 2 also 
maximizes. These rich behaviors can be attributed to 
the high mismatch between the quasi-particle momenta 
on both sides and the interfacial phase difference be¬ 
tween particle and hole channels, which together alter 
the scattering amplitudes in the BTK calculations. 

Finally we discuss the temperature dependence of the 
transport. In Fig. 4 we plot zero-bias NDC vs T/T c ( T c 
is the superconducting transition temperature) at vari¬ 
ous [(a), same convention and setting as Fig. 3(a)] or 
<5y 2 [(b), same as Fig. 3(e)]. The curves can exhibit three 
types of behavior: (1) monotonic increase, (2) monotonic 
decrease, and (3) first increase and then decrease (a reen¬ 
trant phenomenon). These rich behaviors come from the 
same reason as the case of NDC vs V do, because high- 
energy excitations play a more important role at higher 
temperature. As a result, the conductance as a func¬ 
tion of temperature is also sensitive to the bulk SOC and 
spin-orbit-active interface. 


V. CONCLUSION 

In conclusion, a general framework for superconduct¬ 
ing heterostructures with SOC is developed and analyzed 
by generalizing the BTK scheme. Study of the conduc¬ 
tance reveals the crucial role spin-orbit-active interface 
and nonuniform SOC play in determining the transport 
properties. In addition to being directly relevant for un¬ 
derpinning the physics of a range of hybrid systems of 
relevance to various fields including topological matter 
and spintronics, our work opens up many future research 
avenues. One is to apply and extend our transport anal¬ 


ysis to other heterostructures, for example, to study the 
interplay with magnetic fields or ferromagnetism and also 
to study the multichannel interference effects in quasi- 
one-dimensional, higher-dimensional or multiterminal se¬ 
tups [73, 74, 85], also taking into account effects of fluc¬ 
tuations when a low-dimensional superconductor [86] is 
involved. Another is to analyze the impact on pairing 
symmetry, such as the possibility of p -wave pairing [87] 
and its consequences, and to establish a general frame¬ 
work also using complementary approaches such as qua¬ 
siclassic formalism [88]. Finally, another important di¬ 
rection is to explore the consequences of nonuniformity 
and boundary/interface behaviors of SOC in ultracold 
atoms [89] by leveraging their tunability and amenabil¬ 
ity for probing the dynamics. 
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